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Measurement-based quantum computation is a model for quantum information processing utiliz¬ 
ing local measurements on suitably entangled resource states for the implementation of quantum 
gates. A complete characterization for universal resource states is still missing. It has been shown 
that symmetry-protected topological order in one dimension can be exploited for the protection of 
certain quantum gates in measurement-based quantum computation. In this paper we show that 
the two-dimensional plaquette states on arbitrary lattices exhibit nontrivial symmetry-protected 
topological order in terms of symmetry fractionalization and that they are universal resource states 
for quantum computation. Our results of the nontrivial symmetry-protected topological order on 
arbitrary 2D lattices are based on an extension of the recent construction by Chen, Gu, Liu and 
Wen [Phys. Rev. B 87, 155114 (2013)] on the square lattice. 
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I. INTRODUCTION 

Quantum computers provide, for certain classes of 
problems, an exponential speedup over classical comput¬ 
ers [1]. Among many approaches, quantum computa¬ 
tion can be implemented by local measurement on suit¬ 
ably entangled resource states ms, which is known as 
measurement-based quantum computation (MBQC). Re¬ 
source states are said to be universal for MBQC if we can 
implement a universal set of quantum gates by local mea¬ 
surement so that any quantum circuit can be efficiently 
simulated. This model of computation can also be under¬ 
stood in terms of teleportation-based quantum computa¬ 
tion in the valence-bond solid (VBS) formalism [6], as 
well as the so-called correlation-space quantum compu¬ 
tation miH]. A local measurement in MBQC corresponds 
to a joint measurement on partons of the VBS (which are 
to be understood as virtual qubits) within one physical 
site [6] and the resources in the two equivalent models 
are related by a mathematical projection. Since the dis¬ 
covery of the one-way computer [2] on the square-lattice 
cluster state [9], there have been quite a few other uni¬ 
versal resources states uncovered, including, for example, 
graph states HQ], the tricluster state m , weighted graph 
states m, modified toric code states m, and Afffeck- 
Kennedy-Lieb-Tasaki (AKLT) states [13] on various lat¬ 
tices HlHIHI and their deformations mils]. However, as 
of present, there still does not exist a complete charac¬ 
terization of all universal resources. 

An intriguing connection between the resourcefulness 
in MBQC and certain phases of matter was discovered in 
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Ref. [20| , where the authors show that there exists a prop¬ 
erty of many-body states, namely symmetry-protected 
topological (SPT) order in ID, that can be utilized for 
the protection of certain quantum gates in MBQC, even 
though ID quantum states do not accommodate univer¬ 
sal quantum computation. In Ref. [2TI |22| the utility 
of SPT phases for quantum computation in ID has also 
been demonstrated for other symmetry groups beyond 
the Z 2 X Z 2 symmetry. SPT phases are gapped quantum 
systems that have certain symmetries which are not spon¬ 
taneously broken in the nondegenerate ground state [23]- 
l29] . Systems with SPT order (SPTO) can be robust 
against small local perturbations that respect the symme¬ 
try [30l [31] . Ground states of SPT phases can be mapped 
to product states with local unitary (LU) transforma¬ 
tions and thus, are short-range entangled (SRE). Hence, 
they all belong to the same phase as the trivial (prod¬ 
uct) state if no symmetries are imposed. However, if we 
restrict ourselves to symmetric LU transformations, SPT 
ordered systems are distinct from the trivial phase [25^ 
|27|. Due to the potential advantage that SPTO may 
bring to MBQC, one question that arises naturally is 
whether we can find a novel universal resource state that 
exhibits nontrivial SPTO in two spatial dimensions or 
higher. 

In Refs. [28] |29] the authors first discovered a con¬ 
sistent relation between the third group cohomology 
of a symmetry group G and SPTO in 
(2+I)D bosonic systems. Particularly, they prove that 
each nontrivial element of the third group cohomology 
corresponds to a distinct, nontrivial SPT phase. Specif¬ 
ically in Ref. [29], they discuss a fixed point wave func¬ 
tion in its canonical form, showing that it exhibits non¬ 
trivial SPTO with respect to symmetry representations 
constructed from nontrivial elements in PL^iG, U{1)). In 
this approach different SPT phases are not encoded in 
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the wave function but in the different symmetry fraction- 
alizations of the symmetry transformation. The canon¬ 
ical form is defined on a square lattice constituting a 
plaquette-like entanglement structure among partons on 
neighboring sites. In an in-depth discussion, the authors 
of Ref. [29] explain that the canonical fixed point wave- 
function can be related to any other, more generic SRE 
state in the same phase through symmetric LU evolution. 

In this paper, we show that the canonical plaquette-like 
entanglement structure that displays nontrivial SPTO 
enables universal MBQC. Moreover, we shall show that 
such entanglement structure defined on general 2D lat¬ 
tices is a universal resource. But does it also display non¬ 
trivial SPTO as in the square lattice? To answer this, we 
explicitly extend the original symmetry construction to 
arbitrary lattices in two spatial dimensions. With the 
symmetry representation generalized, we can then show 
that these plaquette states and their generalizations on 
all 2D lattices indeed exhibit nontrivial SPTO, depend¬ 
ing on the symmetry. We show that these states are 
universal resources for MBQC as long as their underly¬ 
ing graphs reside in the supercritical phase of percola¬ 
tion. This means that all (generalized) plaquette states 
on 2D regular or quasicrystallene lattices are both non¬ 
trivial SPT states and universal resources for MBQC. 

It is interesting that whether the plaquette states and 
their generalizations are universal resouce states depends 
only on whether the underlying graphs are in the super¬ 
critical phase of percolation. This is exactly the same 
property that governs the universality of graph states, 
which was shown in, e.g., regular lattices m, faulty lat¬ 
tices [32| and random planar graphs [33|. As a contrast, 
the AKLT states, which also exhibit nontrivial SPTO 
only if translation invariance is preserved [29l [34], are 
known to be universal only for some regular lattices m- 

HH]. 


The rest of the paper is organized as follows. In Sec.[n| 
we familiarize the reader with the mathematical frame¬ 
work of group cohomology and its graphic representation. 
Since the results have already been proven in e.g. Ref. 
[29] . we will only outline the most relevant properties. In 
the last part of this section, we will review the aforemen¬ 
tioned canonical form. In Sec. IIP we illustrate the issue 
of defining the symmetry representation on aribitrary 2D 
lattices and provide a solution. We also construct SRE 
plaquette states on arbitrary lattices and verify that they 
exhibit nontrivial SPT order. In Sec. we prove that 
this class of symmetry-protected states is universal for 
MBQC. We summarize and discuss extensions to our re¬ 
sults in Sec. El 


II. REVIEW OF RELEVANT DEFINITIONS 
AND RESULTS 


Symmetry-protected topoplogical phases possess non¬ 
trivial topological order only if a certain symmetry of the 
system is perserved, as opposed to the intrinsic topologi¬ 


cal order for which no symmetry is needed [23H28] . Such 
a SPTO cannot be characterized by a local order param¬ 
eter, either. Ground states of nontrivial SPT phases can¬ 
not be continuously connected to trivial product states 
without closing the gap or breaking the protected symme¬ 
try. From the entanglement viewpoint, these states pos¬ 
sess the so-called short-ranged entanglement as opposed 
to the long-ranged entanglement in the intrinsic topolog¬ 
ically ordered states EH- It turns out that group co¬ 
homology provides a complete classification for ID SPT 
phases [24] and a comprehensive classification for two and 
higher dimensions [29]. Here we review several subjects 
that will be essential for this work, including group coho¬ 
mology and the canonical form for short-ranged entan¬ 
gled states. We follow the seminal paper by Chen, Gu, 
Liu and Wen [29]. We remark that the classification can 
actually be understood via other approaches: (1) anoma¬ 
lous symmetry action at the boundary [35], (2) the non¬ 
linear sigma model with a topological theta term [36l [37] , 
and (3) gauge fields for symmetry twists [38], and (4) 
cobordism [39], all leading to SPT phases beyond coho¬ 
mology [40] . Readers who are familiar with the canonical 
form of short-ranged entangled states in 2D and group co¬ 


for our extension of the original symmetry representa¬ 
tion for SPT phases on arbitrary lattices in two spatial 
dimensions. 


homology can skip this section and go directly to Sec. Ill 


A. SRE many-body states and their canonical form 


In Ref. [29] it was argued intuitively that a SRE state 
in a distinct SPT phase can be tra nsform ed to a sym¬ 
metric tensor network state (see Fig. |l(b)'] ) by symmetric 
LU transformations. It was also argued that this local 
unitary transformation can be accomplished by applying 
a constant-depth quantum circuit. A wave function of 
a generic many-body Hamiltonian is usually very com¬ 
plicated. Fortunately, the above statement allows us to 
formulate a simple SRE tensor network state that can be 
regarded as a fixed-point wave function for more compli¬ 
cated systems in the same phase. Though a generic proof 
is missing. Ref. [29] covers an in-depth analysis and pro¬ 
vides convincing argument, but it will not be repeated 
here in detail. Essentially the argument goes as follows. 

Given a ten sor n etwork representation of any SRE 
state (see Fig. |l(a) ), we can choose squares of size L 
such that entanglement only exists between neighboring 
squares and within. Within each square, we then can 
change the basis and rearrange the degrees of freedom 
using symmetric unitary transformations. As a result, 
we obtai n a t ensor network state in its canonical form 
(see Fig. |l(b)[ ). 

Given a SRE state in its canonical form, we ideally 
want to identify all pairs (1/^, \^)) that exhibit nontriv¬ 
ial SPT order, with Ui being an on-site, global symme¬ 
try representation acting on site i. Therefore, we have 
to clarify how two distinct pairs (I/*, l^p)) and (U[^ W)) 
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FIG. 1. (Color online) (a) Left panel: An exemplary tensor 
network representation of a SRE state. Black dots represent 
degrees of freedom within each site (blue circles). Entangle¬ 
ment exists between connected partons. Together but not 
necessarily individually, black dots form a linear representa¬ 
tion of G on each site, (b) Right panel: Canonical form of a 
SRE tensor network state after rearranging degrees of freedom 
within shaded squares that represent sites in the canonical 
form. 


describing the same phase can be related. There are ba¬ 
sically two equivalence relations: (1) Two phases have 
the same SPT order if their on-site symmetry represen¬ 
tations only differ by a certain class of unitary trans¬ 
formation on each effective site. This is a map 
(Ui^ lip)) ^ (t//, I'lp')) where {ip') has still the same ten¬ 
sor network form. (2) Adding local degrees of freedom 
that form a ID representation of the symmetry does not 
change the phase. The second equivalence relation can 
be utilized to reduce the canonica l form in Fig. |l(b)] to a 
plaquette state as depicted in Fig. 3(a)| Note that it was 
also argued that we can choose to be a symmetry on 
a local region of the canonical form instead of a global 
symmetry. 

One approach to classify SPT order in Ref. [29] is to de¬ 
fine a simple SRE state \i^gs) in its canonical form on pla- 
quettes and find a nontrivial on-site representation of the 
symmetry. Then, more complicated wave functions can 
presumably be found by the aforementioned equivalence 
relations. In this work, we shall mainly focus on plaque¬ 
tte states on arbitrary lattices and construct an on-site 
symmetry representation that is defined on any lattices. 
In Sec.[ni|we will illustrate how the original construction 
of the symmetry operation [29] needs to be modified in 
order to have a symmetry representation defined on all 
lattices and show that this simple modification yields the 
same classification of phases in terms of group cohomol¬ 
ogy* 


B. Group Cohomology 


in Ref. [41]. 

Let ujn • ^ M he a. function of n group elements 

Qi ^ G with value in the G module M. Let C'^{G,M) 
be the space of all such functions which form an Abelian 
group under the function multiplication in M. AG mod¬ 
ule M is an Abelian group with a multiplication opera¬ 
tion compatible with the action of the group G, i.e. 

g • {ah) = {g • a){g -b), \/g eG, a,b e M. (1) 

For our purposes we restrict ourselves to M = U{1) which 
is just the group consisting of complex phases with a 
trivial group action g • a = a for g e G and a G U{1). 

Furthermore, let us define the map dn : C^{G, U{1)) 
C^+i(G,G(l)): 


{dn^n){gi'} 92s •••5 9n-\-l) 

n 

X ■■■,9i-i:9igi+i,gi+2, ■■■,gn+i)- (2) 

i=l 

We utilize this map to define the group of all n- 
coboundaries to be 

B^{G,M) 

— \s^n ^ ^ {G,M)\aJn — 

LUn-l ^ ^(G, M)}, (3) 

and the group of all n-cocycles to be 

Z^{G,M) = {uJn e C^{G,M)\dnUOn = 1}* (4) 

B^{G,M) and Z^{G,M) are Abelian groups satisfying 
B^{G,M) C Z^{G,M). The nth group cohomology is 
the quotient group of n-cocycles and n-coboundaries: 

H^{G, M) = Z^{G, M)IB^{G, M). (5) 

Since we are solely considering two spatial dimensions 
and M = G(l), let us discuss n = 3 as an example. 
From 


{dsUJs){gi,g2,g3s94) 

_ ^3 {92 s93s 94:)(^3 {gi s92g3s 94)^3 {QI s 92 s 9^) 
^3(5^1 ^2, ^3, ^4)^3 (^1 s 92 s 9394 .) 

we obtain the 3-cocyles 

^3 { 92 ,93,9i)^3 {91 > 9293, 54 )^ 3(51 , 92 , 93 ) 


(6) 


= S W 3 


‘^3(5152,53,54)^3(51,52,5354) 


( 7 ) 


As group cohomology shall be used to define the sym¬ 
metry representation, we first give its algebraic definition 
before discussing its graphic representation, with the lat¬ 
ter giving a very useful tool employed in Sec. |III| A more 
rigorous introduction to group cohomology can be found 


and the 3-coboundaries 


bHg,u{i)) 



‘^ 2 ( 52 , 53 )^^ 2 ( 51 , 5253 ) ) 

‘^2(5i52,53)‘^2(51,5253) J 


( 8 ) 
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FIG. 2 . (Color online) Graphic representation of a 3-cochain 
The representation is unique up to a re¬ 
flection along the dotted line. The two mirrored tetrahedrons 
are identihed with opposite chiralities ± 1 . 


the branching structure, we can uniquely determine each 
simplex by assigning opposite chirality ±1 to mirrored 
simplices (see Fig. [^. Given the above two properties, 
we can represent any n-cochain 

( 12 ) 

by a n-simplex with chirality = ±1 and canonically 

ordered vertices gi, gj,g^. 

Even though we can use a n-simplex to represent any 
n-cochain, we have to take into account some redundancy 
implied by the cocycle condition. Since we are predom¬ 
inantly interested in two spatial dimensions, let us only 
consider n = 3 from now on. The 3-cocycle condition 
for cochains takes the form 


which in turn defines the 3rd group cohomology 

nHG, U{1)) = Z\G, U{l))/B\G, U{1)). 

One exemplary representation of a 3-cocycle for a 
group Zd would be (see Ref. M) 

W3(ffi,52,5'3) = exp ^ _ [ 52 + 53 ])^ 

(9) 

for an appropriate c G Z and [g] = g mod d for ^ G 
in the regular representation. 


C. Graphic representation of group cohomology 


92, 93 ,94)1^3(90,gi, 93 ,g 4 )i^ 3 { 9 o,gi,92,gs) _ ^ 

l^sigo, 92,93, 94)’^3{90, §1,92,94) 

(13) 

In particular, it was shown that the cocycle condition im¬ 
plies that two cochains represented by simplices with the 
same surface generate the same phase in M = 1/(1) [29]. 
This statement is true even for a complex with oriented 
loops as long as the simplices within the complex de¬ 
fine appropriate 3-cochains. Due to the one-to-one corre¬ 
spondence between cocycles and cochains, the 3-cochains 
satisfying this condition equivalently define 
On the other hand, 3-coboundaries denoted by A 3 are 
3-cochains that obey 


We now review the graphical representation of the 
group cohomology. It turns out that n-cocycles can 
be represented graphically if we consider so-called n- 
cochains Vn instead. There is a one-to-one relation be¬ 
tween the formerly discussed functions ujn - ^ M 

and cochains. A n-cochain is a map G^~^^ M that 
obeys 


9 ■ ’^n{90,9l, ■■■,9n) = Vn{990,99l, ■■■,99n), (10) 

which is related to ujn by 

(^n{9l,92-,9n) = Vn{l,9l,9l92,-,91 ' ' ' 9n) (H) 

It was shown that any n-cochain can be represented 
uniquely by a branched n-simplex [29]. An arbitrary n- 
simplex can be viewed as a complete graph of its n + 1 
vertices v'^. For example, it is a triangle for n = 2, a 
tetrahedron for n = 3 and a pentachoron for n = 4. Let 
us also impose a naturally ordered labeling of vertices. A 
natural order would imply if i < j. Then, the 

labeling of vertices together with its natural order gives 
rise to a branching structure. A branching structure on 
a simplex is an orientation of edges such that there are 
no oriented loops (see Fig. [^. 

In fact, a branching structure on a simplex uniquely 
determines the canonical ordering of its vertices. On the 
other hand, the branching structure of two mirrored sim¬ 
plices is the same. By imposing a direction in space to 


^ 3 ( 90 , 91 , 92 , 93 ) 


92(91, 92 ,93)92(90, 9 i, 93) 

92(90,92,93)92(90, 9 i, 92) 


(14) 


with /i 2 being a 2 -cochain but not necessarily a cocy¬ 
cle. Then, the equivalence classes [z/ 3 ] := {r' 3 ,z ^3 G 
Z^(G, [/(!))Iz /3 = z^ 3 A 3 } give rise to the third group co¬ 
homology lk^(G,U(l)). 


III. CONSTRUCTING SPT PHASES ON 
ARBITRARY LATTICES IN TWO SPATIAL 
DIMENSIONS 


In this section we first illustrate the issue of defining 
a consistent symmetry representation when generalizing 
the construction in Ref. [29]. It turns out that a simple 
modification can solve this issue and the symmetry can 
be defined on arbitrary lattices, where each site does not 
necessarily have the same number of partons. Moreover, 
the symmetry operation on the boundary encodes useful 
information about the SPTO of the bulk. 


A. The issue of defining a consistent symmetry 
operation 

In accordance with the discussion in Ref. [29], we will 
employ and extend the maximally entangled SRE state 
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\i^gs) used there and examine the global, on-site symme¬ 
try Ui{g) of a d-dimensional finite group G. In fact, one 
key point will be how to properly define such a symme¬ 
try with nontrivial action on the system both without 
and with boundaries. (In fact, examining the bound¬ 
ary will help to understand that the bulk possesses non¬ 
trivial SPTO.) The corresponding exactly solvable, local 
Hamiltonian will just be the projection onto that partic¬ 
ular nondegenerate ground state, in a way similar to the 
CZX model in Ref. [43]. Note that the restriction to fi¬ 
nite groups is useful for the discussion of universality for 
MBQC in Sec. m 

In order to illustrate the issue, let us consider square, 
triangular and honeycomb lattices for now. We define 
the maximally entangled SRE state to be 


I'^gs) — Wpj 

3 

= 0 E •••’ = 9)p, 

3 geG 

(15) 


with parton states |ai), |;d2),..., |C/c) defined on d- 
dimensional Hilbert spaces where group elements label 
physical states according to their regular representation, 
d = \G\ is the order of the group G, Up is the number of 
plaquettes and subscript pj labels the jth plaquette, and 
k G {3,4, 6 } is the number of partons entangled within 
triangular, square and hexagonal plaquettes respectively 
(see Fig. [^. This is a maximally entangled state, i.e., 
the so-called GHZ state, where partons within a plaque¬ 
tte are in the same state in their canonical basis. 

The on-site symmetry on the square-lattice plaquette 
state is defined in Ref. [29] via the 3-cochain z /3 corre¬ 
sponding to Co’s G 1-L^{G^ U{1)) as follows. 


Ui{g)\ai, a 2 , 0 ^ 3 , <^ 4 ) 

= fs{ai,a2,as,a4,g,g)\gai,ga2,gas,ga4), (16) 


where jai, 0 ^ 2 , 0 ^ 3 , 0 L 4 ) is an arbitrary state of partons, g is 
a fixed eleme nt in the group G, (the labeling is illustrated 
in Fig. 4(a)) and 


y*3 (^15 G2 5 0^3 5 ^4 1 91 9^ 

^ a2,9~^g, g)v3{a2, a3,g~^9, S) 

i^3{ai,a4, g-'^g, g)vs{ai, as,g-'^g, g)' 


(17) 


It is shown that such a definition gives rise to a linear 
representation of the group G and the plaquette state is 
indeed invariant under the action of any g on all sites [29] . 
Surprisingly, this construction is easily extended to the 
triangular lattice with six partons per site. In this special 
case, the on-site representation Ui{g) of G G is 

Ui{g)\ai,a2, ..^ae) 

= h{ai,ot 2 , ...,a 6 ,g,g)\gai,ga 2 , -^gaa) (18) 





nZA ( ITl 




nzA ( in 





(a) 





FIG. 3. (Color online) Plaquette states |'0)pj defined on (a) 
square, (b) triangular, (c) honeycomb and (d) kagome lattice. 
Blue dots are identified with partons and shaded circles with 
physical sites. Connections indicate entanglement. 


where jai, 0 ^ 2 ,..., o^e) is an arbitrary state of partons, g is 
a fixed element in G and 


/3 (q!i,Q!2 , ...,ae,g,g) 

^ t^3{ai,a2,g~^9,9)’^3{a2,a3,g~^9,g) 

0 : 4 , g-^g, g)v 3 {oLQ, a^^g-^g, g) 
^ V3(a3,ai,g~^g:g) 
vz{ai.,aQ.,g-'^g,g) 


(19) 


defines the phas e fact or with respect to triangular pla¬ 
quettes (cf. Fig. |4(b)| ). This gives rise to a linear repre¬ 
sentation of G and a global symmetry of analogous 
to the original construction. But how do we generalize 
this to, say, the honeycomb lattice, where each site con¬ 
tains three partons? 


O n any site of one particular sublattice (say, A in 
Fig. 5(a)), a choice of the phase factor is 


fs (^ 15^2 5 ^ 3 :^ 5 ^) 

^ i^3iGi,a2,g~^g,g)ns{a2,a3,g~^g,g) 

us{ai,a3,g-^g,g) 


( 20 ) 


and it will satisfy the group associativity. But on any 
other site of sublattice B if we define the symmetry action 
to be of similar form, then the plaquette state on the 
honeycomb lattice will not be invariant under the global 
on-site symmetry. 

It turns out by careful examination that the choice of 
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FIG. 4. (Color online) (a) Graphic representation of Eq. 

(b) Graphic representation of Eq. ( p^ . Blue dots represent 
partons within one physical site. 


the factor on B sublattice sites as below, 

fs (^1 : ^ 2 : ^35 ^5 

^_ J^3{ai,a2,g~^g,9) _ 

vsiai, a 3 ,g-^g, g)v3{a3, a2,g-^g, g) ’ 


( 21 ) 


will, together with Eq. ( [2Q| ), indeed give rise to a global 
symmetry for the plaquette state. This is because oppo¬ 
site chirality for the graphic representation of a cochain 
is defined by the mirror image and two partons connected 
by a line in Fig. |5(a) are in a maximally entangled state. 
That is, phase factors in Fig. [^arising between parallel 
edges are chosen such that they cancel. We can see that 
whether or not the symmetry representation is a global 
symmetry depends on how we define the on-site branch¬ 
ing structure of fl on sublattice /. In the next section, 
we will generalize this argument to arbitrary 2D lattices. 

Let us point out that just as we can relate the hon¬ 
eycomb lattice to the square lattice by merging vertices, 
we can relate the symmetry representations defined by 
Eq. (20) and (21) on the honeycomb lattice to the rep¬ 


resentation given by Eq. (17) on the sq uare lattice by 
merging sites along edges (^^ Fig. |5(b)'] ), i.e. by forcing 
pairs of merged partons into the same state. 


B. SPT phases from plaquette states on arbitrary 
lattices 


The argument in the previou s sect ion can be extended 
to arbitrary lattices. In Fig. |3(d)| the kagome lattice 
serves as an example for an Archimedean lattice but the 
prescription can be generalized to regular lattices beyond 
this case. In fact, the construction is consis tent o n an ar¬ 
bitrary combination of plaquettes (see Fig. |6(a)[ ) so long 
as (1) neighboring plaquettes share an edge and (2) the 
graph can be embedded on a two-dimensional torus, i.e. 
we impose periodic boundary conditions. 

We start again by defining a maximally entangled state 
\i^gs) on plaquettes around cells on the lattice equivalent 
to Eq. (15) with k being the number of partons within a 



EIG. 5. (Golor online) (a) In order to extend the original 
symmetry representation to the honeycomb lattice, we have 
to define two symmetry representations Uf'{g) and Uf (g) on 
different sublattices A, B. Otherwise, global symmetry is not 
ensured. Red numbers indicate the labeling of partons for 
oi,..., 03 . (b) Merging as a procedure to obtain plaquette 

states on a square lattice from the honeycomb lattice. The 
merging is mathematically a (5-function and is indicated in red 
with respect to the labeling in (a). 




EIG. 6. (Golor online) (a) (Generalized) plaquette states 
\ip)pj on an arbitrary lattice that obeys conditions (1) and 
(2), i.e. that adjacent plaquettes share an edge and that 
we can impose periodic boundary conditions, (b) Perform¬ 
ing measurements in |±)- or Z"^|+)-basis (cf. Eq. 64) on 
virtual qubits that are marked by red crosses, gives rise to a 
valence-bond state. 


plaquette. Note that k is not necessarily constant across 
arbitrary lattices. 

The on-site representation 77/(^) of ^ G G acting on 
site i within the sublattice / is then given by 

e/(fl)|ai, 02, 

= /3(ai,a2, ■■■,ak-,g,g)\gai,ga2, ■■■,9ak-) (22) 

with /c* being the number of partons within one physical 
site, g being a fixed element in G and 

fi{ai,a2, ...,ak-,g,g) 

_ TT ^3{c(ia,C(ia + l,9~^9,9) /23'\ 

where we assume that partons are labeled counterclock¬ 
wise, i.e. i i 1 counting counterclockwise, and 
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FIG. 7. (Color online) Graphic representation of an ex¬ 
emplary phase factor , g,g) as in Eq. For 

the illustrated cutout we find {G} = /c* — 1} and 

{ib} = {2,/c*}. The {aij-plane is illustrated by the shaded 
area. 
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FIG. 8. (Color online) (a) Graphic representation of Eq. 28 
for one choice of {G},{^6}- (b) Graphic representation of 

Eq. I^for the same choice of {G}, {4}- Since the complexes 
have the same surface, they represent the same phase factor 
(cf. Sec. IIC). In fact, for any choice of branching {ia},{4} 
for both complexes, the graphic representations will have the 
same surface. The {a^j-plane is illustrated by the shaded 
area. 


= OLi and The sets of indices 

{^a}, {4} C /c*} define the branching structure of 

/I and are chosen such that {ia} fl {ib} = 0 and U 
{ii)} = {1, On the one hand, {ia} gives counter¬ 

clockwise branched tetrahedrons in the {(a^j-plane of the 
graphic representation of// in Fig. On the other hand, 
{%} defines clockwise orientated tetrahedrons. Natu¬ 
rally, each phase factor (and consequently, each sym¬ 
metry representation) is uniquely specified by the tuple 

In order to illustrate the notation, let us consider the 
kagome lattice as an example. It turns out that we can 
ensure that Uf (g) defines a global symmetry if we divide 
the lat tice i nto three sublattices A^B^C as illustrated 
in Fig. |9(a)[ With respect to the depicted labeling of 
partons, the corresponding phase factors are 

fs (^1 5 •••5 ^ 4 ? ^5 

^ i'3{ai,a2,g~^g,g)v3{a2, « 3 , g) 

vsiai, 0 : 4 , g-^g, g)u3{a4„ as, g~^g, g) 


and 


^ v3{a2,a3-,g~^g,g)v3{a3-,aA,g~^g,g) 

V3 {ol2 , ai,g-^g, g)v3{ai,ai, g-^g, g)' 


(25) 


C. Ui{g) is a linear representation 


In order to show that Ul{g) is a linear representation 
of G, we verify that 

Ul{g'g~^)Ul{g)\ai, a2,afe*) = f// {g')\ai,a2, a**). 

(26) 

Obviously, \aj) ^ W^Lj) \/j for both sides of the equa¬ 
tion. Therefore, let us compare the phase factors. The 
left side of Eq. yields 

TT V3{aia,aia+i,g~^g,g) V3{gai^,gaia+i,gg'~^g,g) 
,, If , V3{ai,+i,ai^,g-^g,g) V3{gai^+i,gai^,gg'-^g,g) 

(27) 

which can be rewritten as 

TT , aia+i,g~^g, g)v3{ai^, ai^+i,g'~^g, g~^g) 

,.-*;f , i^ 3 (a*,+I,a4,5-15,5)2/3(ai.+i,a*,,5'-15,5-15) 

(28) 

since 5 - 2 / 3 ( 50 , 51 , 53 ) = 2 / 3 ( 550 , 551 , 553 ) by definition. 
From the action of Ul{g') we obtain a phase factor 


n 




V3{aig,aig+I,g' ^ 5 , 5 ) 
2/3(a4+i,a4,5'-i5,5) 


(29) 


In agreement with Eq. ( [^ the phase factors can be iden¬ 
tified with tuples (^4, {1)^, {3,4}), (5, {2, 3}, {1,4}) and 
(C, {2, 3}, {1,4}) respectively. In the remaining part of 
this section, we will show that Ul{g) is a linear represen¬ 
tation of G for any {ia}, jib} and that there is indeed a 
choice of I together with a choice of {ia}, {i^} such that 
Ul (g) is a global symmetry of plaquette states defined on 
any lattice. Furthermore, we will see that the symmetry 
action on the system with boundary is the same for this 
model as it is for the original construction. Thus, it gives 
rise to the same classification of SPT phases. 


Comparing the graphic representation of the two phase 
factors reveals that they are the same (see Fig.|^. Hence, 
we have proven that Ul{g) is a linear representation of 
G. This property is completely independent of our choice 
of {ia}, {h} and thus, of the underlying lattice. 

D. Ul{g) is a global symmetry 

To prove that there exists a choice of sublattices I and 
{ia}, {ib} such that (g) is a global symmetry of l^ipgs) 


















(a) 


(b) 


FIG. 9. (Color online) (a) Dividing the kagome lattice into 
sublattices A, B,C as depicted by colored columns and boxes 
together with a labeling of partons as indicated in red, ensures 
that Ui (g) defines a global symmetry if we impose periodic 
boundary conditions between edges along straight lines, (b) A 
unit cell of the kagome lattice. Assuming periodic boundary 
conditions as depicted by dashed lines, this unit cell generates 
the extended lattice by translation. 


on any lattice, we have to show that there exists such a 
choice so that the resulting phase factor F 3 in 

<^Ulig)\^g,)=Fs\^Pgs) (30) 

ij 


equals unity. 

Let us understand why phase factors cancel in the first 
place. Consider two pairs of maximally entangled partons 
as depicted in Fig. We want to choose the symmetry 
representation such that phase factors arising between 
parallel edges cancel. The tetrahedrons representing the 
phase factors have opposite chirality if one can be trans¬ 
formed into the other by reflection. That is, the phase 
factors in Fig. cancel if the branched edges between 
partons are orientated parallelly. As long as we assume 
that all plaquettes share an edge with each neighboring 
plaquette while imposing periodic boundary conditions, 
we can choose a branching between partons such that the 
two graphic representations of phase factors involving all 
four partons can be transformed into each other by reflec¬ 
tion along such edges, i.e. the phase factors cancel (see 
Fig. 10). In other words, there exists a choice of sublat¬ 
tices / and {ia}^{'^b} such that Ul{g) is a global sym¬ 
metry for any lattice that obeys these conditions. This 
is because we can choose w.l.o.g. an arbitrary branching 
on each site. This statement is very general and leaves 
room for a lot of ambiguity, i.e., we may end up choos¬ 
ing different branching structures without respecting the 
periodicity of the lattices. Fortunately, in the case of 
regular lattices, we have to consider only a finite region. 
Therefore, consider the minimal region of some regular 
lattice with plaquettes such that embedding the region 
on a 2 -torus gives rise to the infinite extension of this 
lattice, i.e., the unit cell of a regular lattice. Let us fur¬ 
ther assume that edges of such a region never cut through 
vertices. Such a unit cell contains all structural and sym¬ 
metry information of the macroscopic lattice. In particu- 


(^,g 



I 1^3 ^(7,^,^ ^ g , g ) 


FIG. 10. (Golor online) Partons z+1 and j are in a maximally 
entangled state \f3) and partons i and j + 1 are entangled in 
I7). The graphic representations of exemplary phase factors 
arising between partons is depicted in red. The two opposite 
phase factors cancel if one can be transformed into the other 
by reflection along the edges. This is guaranteed if the orien¬ 
tation of the edges of branched tetrahedrons between partons 
is parallel. 


lar, all edges between vertices of the extended lattice can 
be found within one unit cell. This is the key point since 
the argument for global symmetry as discussed before 
is based on reflections along edges. Given the previous 
considerations, we can find a choice of {^a}, 
on the unit cell such that it defines a global symmetry 
with respect to periodic boundary conditions. By defi¬ 
nition each vertex in the unit cell defines one sublattice. 
Then, {(/, {ia}i {4})} also defines a global symmetry on 
any extended lattice with periodic boundary conditions 
that preserve the lattice configuration. This is because 
any edge configuration of such a macroscopic lattice, can 
be found in the unit cell. That is, an edge in the unit cell 
connecting two vertices of sublattices (/, {ia}, { 4 }) and 
(/', {^ 5 }) respectively, connects the same sublattices 

in the extended lattice in the same way. This is why the 
choice of {(/, {ia}, {h})} ia Eq. (24) and ( [^ gives ri se to 
a global symmetry on the unit cell depicted in Fig. |9(b)| 
and any extension to it as long as it is still embedded 
on a 2 -tor us that preserves the lattice configuration (e.g. 
Fig. |9(a)| ). 

In Fig. we illustrate the branched structures on var¬ 
ious lattices, based on the general proof above. There 
we simplify the drawing and do not show the individ¬ 
ual partons, except in Fig. [If; the arrows represent the 
ordering. 


E. Action of Ul{g) on with boundary 


Now, let us consider the action of llHg) on the sys¬ 
tem \i^gs) with boundary (see Fig. [l^ . We choose 
(/, {ia}i { 4 }) such that the symmetry representation de¬ 
fines a global symmetry. Analogous to the previous anal¬ 
ysis, cochains that depend solely on partons within 
full plaquettes or cochains that involve arguments from 
both full plaquettes and the boundary cancel. However, 
phase factors acting exclusively on the boundary persist. 
Let us relabel the state of plaquettes on the boundary by 
|(ai, 0 ^ 2 ,...) counterclockwise such that \ai)^\aj) are ad¬ 
jacent plaquettes \i \i — j\ = 1. Then, the ground state 
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FIG. 11. (Color online) (a) The unit cell of the triangular 
lattice with plaquettes is indicated in red. Arrows represent 
graphically the on-site action of a symmetry representation 
which is a global symmetry of the corresponding plaquette 
state, (b) Simplified notation for (a), (c)-(h) Exemplary unit 
cells on various lattices together with symmetry representa¬ 
tions that give rise to global symmetries on corresponding 
plaquette states. Each vertex within a unit cell defines a sub¬ 
lattice. 


takes the form 

\^9s) = 

{cxi} 

= Y Y •••)• (31) 


The effective action of (g) on the boundary depends 
on the direction of the bonds between the effective de¬ 
grees of freedom {ai} in Fig. 12 The bonds are iden¬ 
tified with phase factors arising between partons on the 
boundary. The direction reflects the branching of the 
orientated edges in the graphic representation of the 3- 
cochains between such partons. Then, the global effective 
(periodic) action of (g) on the boundary is 


i,I 

= 11*^1" (32) 

{ij) 


where the product runs over all nearest neighbors with 
respect to the direction of the bond and Sij = 1 \i i > j 
and Sij = —1 if i < j. Analogous to the case without 
boundary, we can verify that Ul{g) is a linear represen¬ 
tation on the boundary. 

Eq. (32) coincides with the results in Ref. [29]. That 
is, the action of (g) on the one-dimensional boundary 
can be represented by a matrix-product unitary operator 
(MPUO) U{g), 

u{g)= Y rr(r^’'i(fl)r^’'^(5)---r^’'h5)) 

(33) 

Considering Eq. the matrices T/(^) are found to be 

(Xi+1 

(34) 

for all ai if the bond goes from ai to Other terms 

are zero. If the bond direction is ^ then 

(T]!)"'’®"^^) = Y ’^3{cti+i,cti,9~^9,9)\ai){ai+i\ (35) 

(Xi + 1 


for all ai and other terms are zero. 

Consequently, plaquette states on arbitrary lattices ex¬ 
hibit the same nontrivial properties with respect to the 
symmetry representation Uf (g) as the original model on 
a square lattice. That is, if z /3 is nontrivial, the model 
we constructed has a nontrivial boundary which cannot 
have a SRE, gapped and symmetric ground state. On the 
other hand, if z /3 is trivial, the boundary effective sym¬ 
metry does allow SRE, gapped and symmetric states. 
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(a) 



where qubits are identified with blue dots and sites with 
shaded circles. 

Given the discussion in Sec. Ell we can find an on-site 
representation Ui of the symmetry group Z 2 such that 
this model exhibits nontrivial SPT order. In Ref. [43] it 
is shown that one such representation is 

Uczx = UxUcz (36) 

where 

Ux = Xi 0X2 0X30 X 4 (37) 

Ucz = C Z12C Z23C Z34C Z41 (38) 

wiihCZct = |0)c(0|G)lt + |l)c(l|(8)^t being the controlled- 
phase gate, acting on c(ontrol) and t(arget) qubit and 
X, Z being Pauli X and Z operators. Instead, we could 
consider a symmetry representation constructed from the 
original model of Ref. [29]. Then, Ucz would be substi¬ 
tuted by 


FIG. 12. (Color online) (a) The kagome lattice with exem¬ 
plary boundary. Partons within one plaquette are in the same 
state. The ground state is (b) The boundary 

represented by a ID chain with periodic boundary conditions. 
Bonds represent the action of the symmetry on the effective 
low-energy degrees of freedom |oi, 02 , •••)• 


IV. SYMMETRY-PROTECTED 
TOPOLOGICALLY ORDERED STATES AS 
UNIVERSAL RESOURCES FOR MBQC 

Now that we have shown in the previous section that 
plaquette states on arbitrary lattices are nontrivial SPT 
states, we will now show that these states are also univer¬ 
sal resources for MBQC. We first use the qubit plaquette 
state on the square lattice as an example to illustrate how 
it is a universal resource and then go on to deal with pla¬ 
quette states on arbitrary 2 D lattices. Subsequently, we 
consider the qudit case, which turns out to be a simple 
generalization of the qubit case. Along the way, we also 
discuss how one has to modify the CZX model [43] for 
arbitrary 2D lattices. 


A. 2-level MBQC on the square-lattice plaquette 
state 


Hamiltonian and symmetry. Before we proceed to 
the general case in two spatial dimensions let us consider 
an illustrative example, the CZX model. The CZX model 
was first described in terms of SPT order in Ref. [43] . 
This model has a ground state in the previously discussed 
canonical form and hence, has an entanglement structure 
given by plaquettes, see Eq. (15), represented by squares 


in Fig. 3(a) It is defined on a square lattice with four 
virtual qubits on each physical site as shown in Fig. |3(a)] 


Uscz = sCZi 2 SCZ 23 SCZ 43 sCZi 4 (39) 

where we introduced a special controlled-phase gate 
sCZct = |0)c(0| C + |l)c(l| C If This result assumes 
the specific 3-cocycle representation defined in Eq. © 
with c = 1 and g being the nontrivial element. One can 
easily check that each matrix-product unitary operator 
(see Ref. [29l|43l|44]), representing either one of the two 
symmetry representation on the boundary, corresponds 
to the same nontrivial 3-cocycle. 

The canonical plaquette state is the unique ground 
state of an exactly solvable Hamiltonian H = Hp. 
defined on the jth full plaquette p and the neighb oring 
half plaquettes, i.e. local regions as depicted in Fig. |3(a)| 
The local terms Hp. are given by 

Hp. = -X 4 0 Pj 0 PJ 0 0 Pi (40) 

where X 4 = | 0000 )( 1111 | + | 1111 )( 0000 | acts on a full 
plaquette, P 2 = | 00 )( 00 | + | 11 )( 11 | are projections acting 
on the four surrounding half plaquettes and the Hamil¬ 
tonian acts as identity on the virtual qubits at the cor¬ 
ners. Terms other than X 4 just ensure that Hp. is sym¬ 
metric under the symmetry representations. On the one 
hand, this is because the projections P can effectively 
be treated as 1 with respect to controlled gates not act¬ 
ing on the full plaquette. On the other hand, controlled 
gates acting on P and X4 from the left or the right leave 
the local terms invariant for the same reason for that it 
is a global symmetry. Considering the prescription in 
this paper, this argumentation holds on any other lat¬ 
tice. The Hamiltonian is indeed just the projection on 
the SRE ground state 

\^9s) = ^ (8)(|0000) + |llll))p, (41) 

3 

where rip is the number of plaquettes. 
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FIG. 13. (Color online) (a) An example for two-dimensional 
plaquette states |'0)pj defined in the context of the CZX model 
on a square lattice, (b) Measurement pattern that reduces 
the plaquette state to a bond state. Red crosses mark sites 
on which all virtual qubits are measured in the |di)-basis. 


Quantum computational universality. So far, we 
have predominantly reviewed the construction in [43] . 
Now, let us show that the ground state is indeed uni¬ 
versal for MBQC. In order to process the information 
adequately, we define a local projective measurement pat¬ 
tern that reduces the plaquettes between four qubits on 
a square lattice to bonds between two qubits on another 
square lattice. As we will show, this configuration con¬ 
stitutes a universal resource for MBQC. In Fig. we 
propose a measurement pattern that produces the de¬ 
sired result by measuring each virtual qubit on physical 
qubits along parallel diagonal lines. 

One such strategy to show that the plaquette state 
is a universal resource for MBQC is to construct a local 
measurement that reduces the state to a known universal 
resource state, such as a valence-bond state or a graph 
state. For the square lattice each physical site contains 
four qubit partons and we can express the measurement 
in terms of the degrees of freedom of partons. 

The measurement Mg, used in our first step of reduc¬ 
tion, is a local p rojective measurement on certain sites s 
(see Fig. |13(b)[ ) 


The 1+) outcome projects the other three qubits to (sup¬ 
pressing normalization) 

l(+| • |'04 GHz)i 234 = 1000)234 + 1111)234, (43) 

whereas the |—) outcome projects the other three qubits 
to 


l(—I • |'04 GHz)i 234 = 1000)234 — 1111)234, (44) 


where the latter is the same as the former up to an Pauli 
Z operator on any of the three remaining qubits. By such 
a measurement on all but two qubits, the n-GHZ entan¬ 
glement can be concentrated to any two-qubit standard 
Bell state |00) + |11) (up to possible Pauli correction). 
As illustrated in Fig. measuring sites in a checker¬ 
board pattern in the X basis reduces the entanglement 
structure from plaquettes to bonds between neighboring 
sites given by two-qubit Bell states \B) = |00) + |ll)/v^. 
The resultant state is a valence-bond state (with the bond 
state being \B)) and is, via Eq. (46), equivalent to the one 
defined in Ref. [ 6 ]. The latter was already shown to be 
a universal resource state for MBQC. The above reduc¬ 
tion to a valence-bond state via local measurement can 
be applied more generally to plaquette states on other 
lattices; see, e.g.. Fig. [T4| T he significant difference to 
the measurement in Eq. |42| is that one does not neces¬ 
sarily measure all virtual qubits within one physical site 
(see e.g. Eig. 14(c)[ ). 

To show the universality w e furt her convert the 
valence-bond state shown in Eig. 13(b) to a cluster state 
by local measurement, 

Mc(mi,m 2 ,m 3 ) = H 1 H 2 X (45) 

|nvf'(|oooo)(oooo| + |iiii)(iiii|) IlxfqiJiifs, 


^ i=l 


i=l 


where the Hadamard gate is H = (|+)(0| + |—)(l|)/v^, 
and Hi and H 2 are used to flip the valence-bond state to 
the equivalent one used in Ref. [ 6 ]. In fact, the bonds in 
the two models are related by a Hadamard transforma¬ 
tion: 


4 

= l[Zp\+)i{+\Zl^* (42) 

i=l 

with measurement outcomes rui G {0,1} on virtual 
qubits k = 1 , 2 ,3,4 and |±) = l/\/ 2 (| 0 ) ± | 1 )) being 
the eigenstates of the Pauli X operator. This constitutes 
a valid measurement as = Is is the iden- 

ity operator on site s. The rationale behind this pattern 
is that if a qubit entangled in a n-qubits GHZ state is 
measured in the X-basis, it becomes disentangled while 
leaving the remaining qubits still entangled in a (n — 1 )- 
qubits GHZ state up to an inconsequential local Pauli Z 
byproduct operator in case of an outcome rrik = 1. This 
can be understood with the following example. Eor the 
four-qubit GHZ state |'04ghz) = |0000) + |1111) (where 
we suppress the normalization), a measurement in the X 
basis on the first qubit can have two outcomes: {|±)}. 


H0 1\B) = 10H\B) = (|00) + |01) + |10) - |ll))/2. 

(46) 

The choice of where the Hadamard gates appear is not 
unique, as long as each bond contains one such gate be¬ 
fore the projector. This measurement projects into one of 
the eight two-dimensional subspaces. Regardless of the 
outcome, the resultant state is a cluster state (up to local 
Pauli Z operators), to be proven below. Note that there 
are only three Pauli X operators in Eq. ( [4^ precisely due 
to dividing a 16-dimensional space to 8 two-dimensional 
subspaces. The Pauli X operators can act on any three 
of the four partons. The parameters mi, m 2 , m 3 G {0, 1} 
represent the measurement outcomes. We can define 
the logical zero and one by the mapping (omitting the 
Hadamard gates) 

(|0)(0000| + |1)(1111|)X{^^X^^X^T (47) 
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FIG. 14. (Color online) Illustrated measurement patterns 
that reduce (a) triangular, (b) hexagonal and (c) kagome pla- 
quette states respectively to bond states. Small red crosses 
mark virtual qudits that are measured individually in the |di)- 
or Z"^|+)-basis and big red crosses identify physical sites on 
which all virtual qudits are measured in said basis. 


When mi = m 2 = m 3 = 0 on all sites, the effective 
state in terms of the logical 0 and 1 is exactly the cluster 
state [ 6 ]. We now consider the effect of some mi = 1. 
Due to the commutation relation between Xa and CZab- 

XaCZab = ZkCZabXa, (48) 

an outcome = 1 generates an additional Z byproduct 
operator acting on the qubit at the other end of the bond, 
opposite to parton i\ 


XaCZab\ + +) = Zb CZab\ + +). (49) 

But on any site, if there is an additional Z to the right 
of the projection | 0 )( 0000 | + | 1 )( 1111 |, then 

(|0)(0000| + |1)(1111|)Z, = Z(|0)(0000| + |1)(1111|), 

(50) 

i.e., the effect on the logical qubit can be described by a 
Pauli Z gate. Hence, we can convert the effect of out¬ 
comes mi = 1 to Z gates acting on qubits at the other end 
of the bonds. In other words, the effect is given by logical 
Z gates on the corresponding logical sites. If there is an 
even number of Z gates on a site, the effect is an identity 


{o\ 



(a) (b) 


FIG. 15. (Golor online) (a) Implementation of any one-qubit 
unitary U via joint measurement in basis |iV) on qubits 1,2. 
Qubits 2,3 are entangled in \B). (b) Implementation of the 
controlled-phase gate CZ via joint measurement in basis \0) 
on qubits 1,2,3 and \W) on qubits 5,6,7. Gonnected qubits 
are entangled pairwise in \B). 


gate 1 = Z‘^. On the other hand, if there is an odd num¬ 
ber of Z gates on a site, the effect is a Z gate. Hence, we 
have proven that, regardless of outcomes the post¬ 
measurement state following a measurement (45) on all 


sites is a cluster state up to Pauli Z gates. We remark 
that in the above argument there is an irrelevant overall 
minus sign if there are Pauli X operators acting on both 
ends of the bond. This is due to the commutation of 
S'Hd CZ(2b' 


XaXbCZab = —ZaZbCZabXaXb. 


(51) 


Alternative proof. Another approach to prove uni¬ 
versality after the measurement (42) is to construct a 


universal set of gates using the valence-bond picture as 
in Ref. [ 6 ]. For completeness, we also demonstrate this 
here. Since \B) obviously resembles the bonds in the 
cluster state, the proof of universality is related to the 
argument for universality in the cluster state (see e.g. 
Ref. [6] or Ref. [IS]). That is, we can effectively treat our 
bond state as a cluster state by adjusting the measure¬ 
ment bases. As remarked earlier, the bonds in the two 
models are related by a Hadamard transformation shown 
in Eq. (46). Henceforth, we will omit site labels if it is 


clear from context which site is discussed. 

We already encountered the measurement in the |±)- 
basis that can be implemented to prepare a suitable 
graph from the bond state. Similarly, initialization and 
readout can be performed in the computational basis. 
Application of the theorem in Ref. m gives that any 
one-qubit operation U can be achieved by measuring 
qubits 1 and 2 of a prepared sta te \r])i 0 (| 00 ) + |ll)) 23/2 
in a twisted Bell basis (see Fig. |15(a)| ) 


|V(r,s)) = (C/1'WZ®(8)l)|B). (52) 

with measurement outcomes r, s G {0,1}. This imple¬ 
ments the unitary operation 7/ 17^)3 on an arbitrary input 
state 1?^) = a|0)-}-6|l) up to inconsequential Pauli byprod¬ 
uct operators Z or X in case of outcomes r = 1 and 5 = 1 
respectively. 
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It is well known that we can perform any two-qubit 
transformation by combining one-qubit operations and 
the controlled-phase gate. The latter one can be imple¬ 
mented up to Pauli corrections by performing two dis¬ 
tinct measurements on two adjacent sites in the following 
twisted GHZ-bases (see Fig. |15(b)| ) 

\0{r,s,t)) 

= 1)(W O Z* O X‘)(|000) + |lll))/\/2. 

(53) 

on qubits 1,2,3 with measurement outcomes r, G 
{0,1} and 

\W{u,v,w)) = (X“ (g) Z" (8)X^)(|000) + |lll))/\/2. 

(54) 

on qubits 5,6,7 with measurement outcomes u^v^w G 
{0,1}. It is easy to verify that this does indeed teleport 
and entangle an arbitrary two-qubit input state |x)i ,2 
to (7 Z|x)4^8 up to Pauli operators. This completes the 
alternative proof of universality. 

Note that the logical structure of this proof implies 
that we can perform the quantum computation on the 
plaquette state first and subsequently apply the mea¬ 
surement pattern that reduces the plaquette state to a 
valence-bond state. This procedure works equally well 
since the reductive measurement pattern is independent 
of other measurement outcomes. Hence, such measure¬ 
ments drop out of the temporal order in this model. 


B. 2-level plaquette states on arbitrary 2D lattices 


Here we show that, in addition to the square lattice, 
plaquette states on arbitrary 2D lattices are indeed uni¬ 
versal resource states for MBQC. By lattices we mean 
arbitrary 2D graphs that are polygon triangulation of 
two-tori with a macroscopic number of vertices and edges 
(as well as faces). (Note polygons here need not be con¬ 
vex, as long as they are topologically equivalent.) This 
definition also implies that the graphs we consider for 
universality are in the supercritical phase of 2D perco¬ 
lation. That is, all the regular 2D lattices, such as the 
eleven Archimedean lattices and their dual lattices, host 
plaquette states that both exhibit nontrivial SPT order 
and serve as universal resource for quantum computation. 


Hamiltonian and symmetry. We can also generalize 
the CZX model to other lattices. Hamiltonians similar 
to the one introduced in Eq. (40) can be defined on ar¬ 
bitrary lattices so that the unique ground states are pla¬ 
quette states on the same lattices. However, if the vertex 
degree is odd (such as in the honeycomb lattice) then the 
symmetry generated by Uczx = UczUx is no longer a 
linear representation of Z 2 , as for example. 


Ux = X 1 X 2 X 3 , Ucz = CZ 12 CZ 23 CZ 31 , (55) 

^czx = (56) 


due to UczUx = —UxUcz when there is an odd 
number of partons. But as remarked in a footnote of 
Ref. [43] this can be remedied by replacing X by iX^ 
i.e., using Uix = {iXi){iX 2 ){iX 2 ,) and the iCZX group 
{I, Uiczx = UczUix] gives a Z 2 symmetry group of the 
Hamiltonian. Alternatively, the special CZX unitary can 
defined consistently and is a symmetry of the Hamilto¬ 
nian. This is because it is guaranteed that 

UsCzUx = UxUscz (57) 

since X^Xt sCZ^tX^Xt = ZcZtsCZct instead of 
X,XtCZ,tX,Xt = -Z,ZtCZ,t. 

Quantum computational universality. In the previ¬ 
ous section we discuss the concentration of GHZ entan¬ 
glement in a plaquette to a Bell state between any two 
partons on this plaquette. Schematically, we can use a 
dashed line connecting these two partons (or just the two 
sites containing these two partons) to represent the said 
concentration. This represents an operation (which we 
shall refer to as 01) that can be implemented by local 
measurement; see Fig. [l^a) for illustration. 

Moreover, suppose we have two pairs of 4-GHZ states, 
|'04GHz)i,2,3,4 and |'04GHz)a,6,c,d and the qubits I and a 
belong to partons on the same site. Then, it is worth 
mentioning here that a measurement of I and a in the 
Bell-state basis {|00) ± |II), |0I) ± |I0)} (which is still 
local on the site containing I and a) can merge the two 
GHZ states into one shared among qubits 2, 3,4, 6, c, d, 
for example. 


la((00| ± (II|) • |'04GHz)l234'04GHz)a6cd 

= |000000)2346cd ± |IIIIII)2346cd. (58) 

Other outcomes |0I) ± |I0) give equivalent merging up to 
Pauli Z and X operators. 

Suppose further that qubits 2 and b also belong to one 
other site that shares two plaquette edges with the site 
containing qubits I and a. Then, as illustrated schemat¬ 
ically in Fig. [l^b), the two plaquettes are merged into 
one. (For completeness, note that one can also perform 
a measurement on qubits 2 and b to remove both qubits 
from the plaquette or just on one of them to keep the 
other parton.) Together with the previous concentration 
operation, the entanglement of two plaquettes (or more) 
can be concentrated to one Bell state between any two 
parton qubits on any two separate locations of the two (or 
more) plaquettes. One can schematically draw a dashed 
line connecting these two partons or, for simplicity, the 
two physical sites containing the two partons, as shown 
in Fig. [T^b) and later in Figs. 17 and 18 This repre¬ 
sents another operation, which we shall refer to as 02, 
that can be implemented by local measurement. 

The above schematic operations 01 and 02 are suffi¬ 
cient to show that plaquette states on arbitrary 2D lat¬ 
tices are universal resources for MBQG. (In fact 01 is 
already sufficient, but the inclusion of 02 will be use¬ 
ful.) If the entanglement of a plaquette is concentrated 
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FIG. 16. (Color online) Illustration of operations (a) 01 and 
(b) 02 as described in the text. For the 01 operation, entan¬ 
glement can be concentrated to any two partons on the pla- 
quette. For the 02 operation, two plaquettes can be merged 
into one single plaquette and entanglement can be concen¬ 
trated to any two partons on these two plaquettes. 



to two partons, no other partons can be entangled. This 
implies that there cannot be a second dashed line drawn 
inside the same plaquette or polygon (i.e., there can be 
at most one), nor can there be more than one dashed 
line across any adjacent plaquettes. The key point of our 
proof is that among arbitrary 2D graphs in the super¬ 
critical phase, one can always find a minor graph that 
is topologically equivalent to a honeycomb lattice [32] . 
with vertices separated from one another by a few graph 
distances and their number being only a fraction of the 
original number of vertices. If they are far enough apart, 
one can always connect these vertices with dashed lines 
cutting through polygons such that there is at most one 
line inside each polygon. This is illustrated in Fig. 
Thereby, we obtain a valence-bond state on a honeycomb 
lattice, which can be further converted to a cluster state 
on the same lattice by local measurement, simply ex¬ 
tending the proof described earlier in Sec. |IV A| for the 
square-lattice valence-bond state to the honeycomb case. 
Thus, we have proven that plaquette states on arbitrary 
2D lattices (those whose graphs are in the supercritical 
phase of percolation) are universal resources for MBQC. 

As a remark, partons within plaquettes that do not 
participate in the entanglement concentration (i.e., in 
those plaquettes where no dashed lines were drawn), can 
be measured in any basis so as to remove them from 
the resource state and computation. We illustrate the 
reduction procedure from plaquette states to valence- 
bond states on various Archimedean lattices in Fig. 

As one can see, the operation 01 is sufficient. But for 
convenience, dashed lines corresponding to an operation 
02, i.e. those that cross two polygons, were used in 
Fig. ^d). They could be replaced by combinations of 
01 dashed lines. The depicted prescriptions are not nec¬ 
essarily optimal but are used to demonstrate the proof 
described above. 


FIG. 17. (Golor online) Schematic illustration for the proof 
that plaquette states on arbitrary lattices are universal re¬ 
sources. Solid dark blue curves represent paths traversing 
edges on the lattice and dashed red curves represent concen¬ 
tration of entanglement via faces or polygons (which we did 
not draw explicitly). If the vertices where blue curves cross 
are sufficiently far apart, it is guaranteed that no two red 
dashed lines are to cross the same face (or polygon). Hence, a 
valence-bond state on a honeycomb lattice can be generated 
via local measurement. 


C. d-level MBQC on plaquette states 


Originally, cluster states were defined for 2-level sys¬ 
tems [9] but they have been generalized to d-level systems 
for MBQC as well [45l|47]. Here, we generalize the results 
of the previous sections for qubit plaquette states to d- 
level systems with ground states \'ipgs) given by Eq. (15), 
i.e. with virtual qadits replacing qubits: 


I'^gs) — 


1 


- ^ n ,/2 



d-1 


=9,132 = g,...,Ck = g) 


j \^=0 


Pj ’ 


(59) 


which is constructed from the prescription in Sec. |m|for 
a finite d-dimensional symmetry group on any plaquette, 
i.e. arbitrary k. 

Let us first consider the square lattice, that is, k = 
k^ _ 4 measurement M applied in this case is just 
the extension of Eq. (42) from qubit to qudit systems. 


k* 

= (60) 

i=l 

where Z^’s are diagonal matrices (generalizing the Pauli 
Z operator) acting on the ith virtual qudit within site s 

Zi = diag(l, Q, ..., Q^“^) (61) 
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FIG. 18. (Color online) Various lattices for the reduction of 
plaquette states to valence-bond states. Utilizing the method 
described in the text and illustrated in Fig. |16| plaquette 
states on black lattices are reduced to valence-bond states 
on lattices depicted by dashed red lines. 


with Q = exp(^) such that = 1. 


l+)i 


Vd ^ 

V « j=0 


d-1 


(62) 


is an eigenvector of the generalized Pauli X operator X 


that acts as 


Xi\j)i = \j + 1 mod d)i. 


(63) 


X and Z generate the generalized Pauli group or the 
Heisenberg-Weyl group on a single qudit system. To 
show that the generalized plaquette states are also uni¬ 
versal resources, we will construct a similar measurement 
that reduces them to qudit valence-bond states. The 
valence-bond picture of MBQC [6] was generalized to qu- 
dits in Ref. [45]. The MBQC with qudit cluster states was 
also discussed in m- As for the qubit case universality 
can be proven in two different approaches. 

It is easily verified that Eq. (60) is a complete measure¬ 


ment that disentangles the ground state in Eq. (59) anal¬ 


ogous to the measurement in the preceding section. The 
resulting state is a valence-bond state with the two-qudit 
bonds being \B) = bi)/v^ up to inconsequential 

Pauli operates. 

Analogous to the discussion in Sec. |IVB 


we can ex¬ 


tend the measurement (60) to arbitrary 2D lattices. Eor 


some simple lattices, we can find a similar measurement 
pattern that reduces systems defined on plaquettes to 
valence-bond states. In doing so, we can reduce, for ex¬ 
ample, the plaquette state on a triangular lattice to a 
bond state on a hexagonal lattice or the plaquette state 
on a honeycomb lattice to a bond state with the same 
hexagonal periodicity (see Eig. 14). 

Let us point out that such a measurement pattern is 
chosen due to its simplicity but one can think of pat¬ 
terns that reduce the square lattice to more complicated 
structures such as e.g. a decorated square lattice (which 
includes additional sites on edges of the square lattice). 
In that case, we do not necessarily measure each virtual 
qubit within one physical site and the post-measurement 
state retains some randomness on such sites. That is, the 
measurement Mg on site s acts on virtual qubits A, •••Gz 
for some I < /c* according to 


k=l 

with outcomes rrii-^^ G {0, ...,(i — 1} and acts as 

unity everywhere else. Such a measurement pattern is 
necessary in order t o reduc e most lattices to valence-bond 
states (see e.g. Eig. 14(c)[ ). It has already been employed 
with d = 2 in the context of operations 01 and 02 in the 
previous section. In fact, the schematic diagrams for this 
reduction are the same, regardless whether partons are 
qubits or qudits; see Eig. 16 Considering Eq. (64), it 


is intuitively clear that we can reduce any generalized 
plaquette state to a general ized v alence-bond state by 
local measurement (see Eig. |6(b)[ ). The schematic dia¬ 
gram for this proof is the same as for the qubit case; see 
Eig. Hence, generalized plaquette states defined on 
arbitrary lattices can be reduced to generalized valence- 
bond states. 

Analogous to the qubit case, we can employ local mea¬ 
surement to convert generalized valence-bond states to 
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cluster states. For example, the generalized valence-bond 
state on the square lattice can be reduced to a cluster 
state via local measurement, similar to Eq. 

Md{mi,m 2 ,ms) = x (65) 

I n (E n (^0"' I F1F2, 

e=i ■J'=o *=i J 

where the Hadamard gate of the qubit case is substituted 
by a quantum Fourier transformation F: 

d-l 

(66) 

k=0 


We note that, similar to the qubit case, the choice of 
where the F gates appear is not unique, as long as each 
bond contains one such gate before the projector; see also 
Eq. (72) below. This measurement projects into one of 
the crd-dimensional subspaces. Regardless of the out¬ 
come, the resultant state is a d-level cluster state (up to 
local Z operators), to be proven below. The parameters 
mi, m 2 , m 3 G ,d— 1 } represent the measurement 

outcomes. Analogous to the qubit case, we can define the 
logical state by a mapping (omitting the F gates) 


(E \j){mmrxpxr^. (67) 

j=0 

When mi = m 2 = m 3 = 0 on all sites, the effective 
state in terms of the logical j is exactly the d-level cluster 
state m- We now consider the effect of some m^ 7 ^ 0. 
The commutation relation between the d-level Xa and 
CZab is 


XaCZab = Z-^CZabXa, (68) 

where the d-dimensional CZ-gate acting on c(ontrol) and 
t(arget) qudit is given by 

CZ\jk)c,t = ^^''\jk)c,f (69) 

From the above we see that an outcome m^ 7 ^ 0 generates 
an additional byproduct operator acting on the 

qudit at the other end of the bond, opposite to parton i: 

17“ CXb|++) = CZab\++). (70) 


gates on the corresponding logical sites. Hence, we have 
proven that, regardless of the outcomes m^, the post¬ 
measurement state following a measurement ( [4^ on all 
sites is a d-level cluster state up to Z gates (to some 
power depending on measurement outcomes). 

Alternatively, we can also show that the post¬ 
measurement state subsequent to a measurement (64) 


is universal for MBQC by explicitly constructing a uni¬ 
versal set of quantum gates by local measurements. In 
Ref. [ 45 ] it has been shown that a d-level cluster state is 
universal for MBQC. Similar to the qubit case, the bonds 
\B) and the bonds in the cluster state are related by 


d-l 


{F^1)\B) = (1. 


= E (72) 

^ ^ k,j=0 


Note that we ignored and will continue to ignore site 
labels. Due to the relation in Eq. (72), we can again twist 
the measurement bases and projections by F in order to 
show that the valence-bond state is universal for MBQC. 

The measurement in the Z’^|+)-basis can be deployed 
for preparation. Once again, measurements in the com¬ 
putational basis serve for initialization and readout. 
Since the proof in Ref. m is valid for any dimensionality 
d < 00 of a single parton’s Hilbert space, we can easily 
extend the measurement basis in Eq. [^to qudits: 

d-l 

\N{r, s)) = (U^X^Z^ 0 1)(^ \jj))/Vd (73) 

j=0 

with measurement outcomes r, s G {0,1,..., d — 1}. Uni¬ 
versality is ensured if we can combine any one-qudit op¬ 
eration with the CZ-gate (see Ref. [45]). The CZ-gate 
is implemented up to Pauli byproduct operators via two 
distinct measurements on two adjacent sites in the bases 


\0(r,s,t)) 


d-l 


= (1 0 Ft 0 1 )(X’- 0 0 * 0 X*)(^ \jjj))/Vd (74) 

7=0 

on qudits 1,2,3 with measurement results r, G 
{0,1,..., d — 1} and 

d-l 

\W{u, V, w)) = (X“ 0 0 X"')(^ \jjj))/Vd (75) 

7=0 


But on any site, if there is an additional Z to the 
right of the projection \j){jjjj\^ then 

(E (71) 

7=0 7=0 

i.e., the effect on the logical qudit can be described by 
a Z gate. Hence, we can convert the effect of outcomes 
7 ^ 0 to Zi gates acting on qudits at the other end of 
the bonds. In other words, the effect is given by logical Z 


on qudits 5,6,7 with measurement results u^v^w G 

{ 0 , 1 ,-V-i}. 

This completes the proof of universality for the general 
d-dimensional case. Since the argument is independent 
of the underlying lattice, we can utilize any generalized 
plaquette state for quantum computation. 

Note that the logical structure of our proof implies 
again that we can perform the quantum computation on 
the plaquette state first and subsequently apply the dis¬ 
cussed measurement pattern M accordingly. This pro¬ 
cedure works equally well since the measurement M is 
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independent of other measurement outcomes. Hence, 
such measurements drop out of the temporal order in 
this model. 


V. CONCLUSION 

In this paper we extend the symmetry representation 
of the group cohomology approach for SPTO by Chen 
et al. [29] to arbitrary lattices in order to show that 
the plaquette states exhibit nontrivial SPTO not only 
on the square lattice but also on arbitrary 2D lattices. 
In connection to quantum computation, we prove that 
plaquette states and their generalizations serve as uni¬ 
versal resource states for measurement-based quantum 
computation as long as the underlying graphs reside in 
the supercritical phase of percolation. It is interesting to 
note that the deciding property for universality in graph 
states is exactly the same [T0l[32l|33]. As a contrast, the 
AKLT states, which also exhibit nontrivial SPTO only if 
translation invariance is preserved [29l [34] , are known to 
be universal only for some regular lattices HMH). 

Our initial motivation was to see if one could reach 
similar conclusions as in the ID case by Else et al. m, 
i.e. that certain quantum gates (or even better, a set of 
universal gates) may be protected within the entire 2D 
SPT phase (at least for a certain symmetry). If this was 
the case, then SPT phases would provide naturally pro¬ 
tected resource states and could potentially yield better 
fault tolerance in comparison to other, topologically triv¬ 
ial universal resources. Measurements would then only 
needed to be performed in certain (symmetric) bases, 
giving rise to the protected gates and possibly reducing 
the types of error. But we have not succeeded in doing 
so. As for the plaquette states, the measurement that 
reduces them to cluster states can be done at one step, 
so the fault tolerance threshold the same as that of 2D 
cluster states. Unfortunately, this is of order 10“^ to 
10“^. But for 3D cluster states, the threshold is as high 
as 10“^, since the measurement patterns can simulate 
braiding of anyons and the topological protection can be 
exploited [48| . Similarly, one can use a three-dimensional 
version of plaquette states and reduce them to 3D clus¬ 
ter states. This means that one can achieve the same 
threshold as for the 3D cluster state. However, we do not 
know whether quantum computation performed directly 
on plaquette states without reduction to the cluster state 
could result in a better threshold. 

Similar to 2D graph states, the plaquette states can 


be defined on any 2D graph. However, the local Hilbert- 
space dimension can be large, depending on the number 
of partons on each site. Thus, it would be desirable to 
find projections from partons to a smaller physical sub¬ 
space while at the same time preserving the nontrivial 
SPTO and quantum computational universality. More¬ 
over, we have solely demonstrated the universality of pla¬ 
quette states, which are nothing but representative states 
in SPT phases. One immediate question is whether all 
ground states in any or some nontrivial SPT phase are 
still universal. These are open-ended questions for future 
investigation. 

In this paper, we solely discuss bosonic systems. How¬ 
ever, it has been shown in Ref. [43] that SPT order exists 
for fermionic systems in 2D as well. In fact, plaquette 
states play a similar role if fermions are under consid¬ 
eration. In Ref. [49] the authors argue that SPT order 
in some fermionic systems can be classified by a special 
group supercohomology theory. Furthermore, it has al¬ 
ready been shown in Ref. m that fermionic MBQC is 
possible provided that we consider additional information 
processing channels, namely parity modes, besides the 
conventional information mode. The additional modes 
carry the parity property of fermions. We suspect that 
a fermio nic sys tem defined on e.g. triangular plaque- 
ttes. Fig. |I4(a)[ exhibits nontrivial SPT order and that 
it might also serve as a fermionic resource state. 

As we were finishing up this work, we became aware of 
a related work by Miller and Miyake in Ref. m where 
the authors construct a specific symmetry-protected 
qubit state and demonstrate that it is universal for 
MBQC. 
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